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ABSTRACT 

Experimental  and  analytical  studies  are  made  of  the  elastic  and 
elastic-plastic  response  of  a  three  degree  of  freedom  dynamical  system 
composed  of  three  rigid  and  compact  masses  connected  by  slender  columns, 
the  bottom  (heaviest)  mass  being  suspended  as  a  ballistic  pendulum  and 
subjected  to  an  initial  velocity  step.   Such  a  system  is  related  to  the 
hehavior  of  machinery  mounted  in  a  ship  subjected  to  an  underwater  ex- 
plosion. 

The  analytical  work  was  performed  on  a  high  speed  digital  computer 
and  employed  mathematical  description  of  the  structural  members  in  which 
force-deflection  relations  depended  upon  previous  plastic  strain  history. 

Comparisons  are  made  between  experimental  and  analytical  responses 
on  the  basis  of  the  velocity-time  history  of  the  two  smaller  masses. 
Good  agreement  is  shown. 
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duct  ion 

The  present  interest  in  the  response  of  structures  subjected  to 

rapn'  lly   pplied  dynamic  loads  has  produced  an  increasing  amount  of 

literature  in  the  area  of  ly      loading.   The  literature  covers  the 

behavior  of  structural  metals,  the  res]  pie  structures,  and  the 

response  of  ships  subjected  to  ex- ]        111  ,  I  2 ■   , 

A  recenl    /estigation  c  "he  U. 

Naval  Fostgraduate  Schoo      I  erey,  California,  dealt  with  the  behavior 
M  pie  plane  structure  -     L  of 

great  stiffness  compare."  to  the  col  Lch  the  greatest  part 

of  the  mass  and  weight  was  i  The  loa       s  an  impact  in 

lane  of  the  structure  and  at  the  roof  level.   This  system  could  be 
adequately  approximated  by  a  stngle  degree  of  freedom  model  and  a  simple 
theoretical  analysis,  based  essentially  on  energy  considerations, 
able  to  make  good  predictions  of  experimental  behavior. 

In  the  case  of  a  similar  multi-story  structure,  the  energy  is  divided 
between  the  various  roof  an    '     asses  and  it  is  not  possible  to  use 
energy  considerations  alone  to  obtain  a  theoretical  prediction  of  dynamic 
behavior.   Thus,  analysis  of  a  system  abo      \ngle  degree  of  freedom 
requires  evaluation  of  the     '       ations  of  motion  of  the  system.   If 
the  structure  passes  into  tl      stic  region,  a  suitable  assumption  must 
be  made  to  compute  the  behavior  in  this  region.   The  numerical  manipula- 
tions required  in  solving  such  a  problem  become  enormous  and  the  use  of  a 
high  speed  digital  computer  becomes  necessary. 

This  thesis  deals  with  a  small  two-story  structure  mounted  on  a 
large  movable  mass.   The  entire  system  was  reduced  to  a  three  degree  of 


^Numbers  in |  ]  refer  to  Bibliography, 


freedom  system.   This  particular  system  can  be  likened  to  a  piece  of 
equipment  mounted  to  a  ship's  hull.   The  ship  is  then  subjected  to  an 
underwater  explosion  and  the  system's  response  is  studied.   Reference  I  3  I 
gives  examples  of  equipment  and  machinery  damage  and  response,  wherein 
illustrations  of  actual  shipboard  foundation  failures  are  presented  along 
with  a  typical  velocity  history  for  a  mild  shock.  The  ability  to  pre- 
dict the  response  of  such  a  system  would  aid  in  the  original  design  of 
such  structures. 

A  mathematical  model  was  evolved  for  the  above  system  and  a  compari- 
son made  between  the  response  predicted  by  the  mathematical  model  and  a 
test  structure. 

It  is  felt  necessary  to  mention  the  fact  that  the  experimental  por- 
tion of  the  project  was  carried  out  prior  to  the  evaluation  of  the  mathe- 
matical model  and  certain  complications  arose  in  the  mathematical  evolu- 
tion that  would  have  been  eliminated  had  the  sequence  been  reversed.   The 
values  obtained  from  the  mathematical  model  will  be  referred  to  as  the 
"theoretical  results"  throughout;  however,  these  results  cannot  be 
truly  classified  as  "theoretical"  since  a  large  amount  of  experimentation 
with  the  mathematical  model  was  necessary  before  a  suitable  comparison 
could  be  made. 

The  experimental  portion  of  the  thesis  was  conducted  entirely  at 
the  David  Taylor  Model  Basin,  Carderock,  Md.,  during  the  period  30  June  to 
19  July  1961.   The  theoretical  investigation  and  reduction  of  data  was 
carried  out  at  the  U.  S.  Naval  Postgraduate  School,  Monterey,  California, 
during  the  period  1  January  1962  to  1  April  1962,  utilizing  the  Control 
Data  Corporation  1604  computer. 


The  experimental  investigation  will  be  discussed  before  the  theo- 
retical considerations  due  to  the  fact  that  this  was  the  order  in  which 
they  were  actually  conducted.   A  better  understanding  of  the  complications 
introduced  in  the  theoretical  considerations  will  be  obtained  if  the  factors 
responsible  for  the  complications  are  presented  first. 


2.   Experimental  Investigation 
Test  Structure 

The  test  structure  consisted  of  two  rectangular  hot  rolled  steel  bars 
measuring  2"  x  2"  x  12",  separated  and  supported  by  two  sets  of  1/8"  x  2" 
hot  rolled  steel  columns.   Figure  (la)  shows  a  sketch  of  th  structure  with- 
out the  attached  measuring  equipment.   The  only  variation  made  to  the  test 
structure  was  the  introduction  of  3003-H14   aluminum  columns  for  ten  (10) 
of  the  twenty-six  (26)  test  runs.   These  columns  also  measured  1/8"  x  2". 

The  length  of  the  supporting  columns,  prior  to  mounting  was  12". 
Holes  measuring  3/8"  in  diameter  were  drilled  on  the  centerline  one  inch 
from  each  end  of  the  supporting  columns  to  accommodate  a  3/8"-NC  bolt  for 
clamping  the  supports  to  the  base  and  bars  mentioned  above. 

The  rectangular  bars  and  base  were  drilled  and  tapped  for  3/8"-NC 
threads  3/4"  deep  at  both  ends,  to  receive  the  securing  bolts.   Holes 
measuring  3/8"  in  diameter  were  also  drilled  in  the  sides  of  the  bars  for 
attachment  of  the  velocity  meters  used  in  the  test.   Small  metal  bases, 
fitted  to  receive  the  accelerometers  used  in  the  test,  were  welded  to  the 
upper  surfaces  of  each  of  the  rectangular  bars. 

The  entire  clamping  of  the  supporting  columns  was  accomplished  by 
the  single  3/8"-NC  bolt  along  with  a  1/8"  x  2"  x  2"  clamping  plate,  see 
figure  (lb).   The  spacing  between  the  base  to  the  center  bar  and  the 
distance  between  the  center  bar  to  the  upper  bar  was  approximately  8". 

The  single  screw  attachment  of  the  columns  to  the  base  and  bars 
proved  to  be  weak  as  far  as  giving  good  fixity  at  these  junctions.   This 
poor  fixity  caused  the  effective  length  of  the  columns  to  be  somewhat  in 
excess  of  the  actual  free  length  of  8".   Figure  (2)  is  a  picture  of  the 


fully  instrumented  test  structure.   The  velocity  meter  coils  are  on 
the  side  seen  on  the  picture.   The  accelerometers  are  on  the  upper  sur- 
faces facing  toward  the  left. 

The  base  of  the  structure  contained  two  2"  x  2"  x  2"  hot  rolled  steel 
blocks  to  be  used  for  attachment  of  the  supporting  columns.   The  base  was 
designed  to  fit  a  mounting  block  already  attached  to  the  ballistic  pendu- 
lum used  in  the  tests. 

Since  the  width  of  the  supporting  columns  was  large  in  comparison  with 
the  thickness,  the  structure's  motion  was  limited  to  a  single  plane. 

Ballistic  Pendulum 

The  ballistic  pendulum  was"  used  as  the  means  for  producing  the 
impulsive  loading  on  the  structure.   A  picture  of  the  pendulum  is  shown  in 
figure  (3). 

The  ballistic  pendulum  consisted  of  four  major  components;   1.  The 
anvil,  upon  whicb  the  test  structure  was  mounted.   2.   The  hammer,  which 
transferred  its  energy  to  the  anvil  producing  the  impulsive  load.   3.   The 
hammer  hold-back  and  release  mechanism.   4.   The  anvil  hold-back  and 
catching  system. 

The  hammer  consisted  of  a  solid  steel  cylinder^  8"  in  diameter  and  48" 
in  length.   Four  smaller  cylinders  were  mounted  on  the  sides  to  receive 
the  supporting  wires  from  the  overhead.   These  cylinders  were  4"  in  d  a  eter 
and  6-1/2"  in  length.   The  total  x^eight  of  the  hammer  was  785.1  lbs. 

The  anvil  was  likewise  8"  in  diameter;  however,  its  length  was  36". 
The  four  supporting  cylinders  measured  3.5"  in  diameter  and  7"  in  length. 
The  supporting  cables  were  8  ft  in  length  ±1".   The  total  weight  of  the 
anvil  was  598.6  lbs. 
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d  release  mechanism  consisted  of  an  electric 
hoist  I     11  the  hanui  rated  position,  the  initial  posi- 

'.'.pendii,  rhe  rel  'nism  consisted 

of  an  electro-magnetic  latch  which  could  be  operated  from  a  remote  station 
near  the  recording  instruments. 

The  anvil  hold-back  and  catching  system  cor.sis.ted  of  a  ratchet  pulley 
and  a  tie-back  line.   The  line  was  attached  to  the  anvil  at  one  end,  passed 
through  the  pulley  and  carried  a   25  lb,  weight  at  the  free  end.   Upon  im- 
pact the  anvil  was  forced  back  and  elevated.   The  25  lb4  weight  pulled  the 
line  through  the  ratchet  pulley  and  the  anvil  was  prevented  from  return- 
ing upon  reaching  its  maximum  elevation. 

The  ratchet  catching  system  was  used  to  prevent  interference  between 
the  anvil  supporting  wires  and  the  test  structure.   The  proximity  of  the 
test  structure  and  the  anvil  supporting  wires,  when  the  anvil  was  in  its 
equilibrium  position,  was  such  that  the  structure  came  into  contact  with 
the  supporting  wires  upon  impact.   To  eliminate  this,  the  anvil  was  pulled 
back  approximately  8-1/2"  horizontally  prior  to  impact. 

Since  the  experimental  testing  was  conducted  prior  to  a  thorough  theo- 
retical investigation,  the  significance  of  this  initial  anvil  position  was 
not  realized,  and  the  exact  initial  horizontal  location  of  the  anvil  was 
not  recorded.   Upon  examination  of  the  theoretical  aspects  it  was  found  that 
this  initial  position  was  quite  significant  and  complicated  the  theoretical 
solution  of  the  problem.   More  will  be  said  concerning  this  later. 

Measuring  Equipment 

The  test  structure  was  instrumented  with  three  bar-magnet  velocity 
meters  manufactured  and  calibrated  by  the  D-^vid  Tnylor  Model  Basin.   One 


velocity  meter  was  mounted  on  each  of  the  three  concentrated  masses  of 
the  structure.   The  meter  mounted  on  the  anvil  had  a  sensitivity  of  146 
mv/fps,  the  meter  on  the  center  mass  (lower"  rectangular  bar)  had  a  sensi- 
tivity of  167  mv/fps,  and  the  third  meter,  mounted  on  the  upper  mass  (upper 
rectangular  bar),  had  a  sensitivity  of  151  mv/fps.   From  tests  conducted  by 
the  David  Taylor  Model  Basin,  it  was  determined  that  the  sensitivities 
varied  approximately  ±3%  over  a  4  inch  stroke. 

Two  Statham  accelerometers  were  mounted,  one  each,  on  the  upper  and 
center  masses  of  the  test  structure. 

The  output  from  the  velocity  meters  went  to  the  Velocity  Meter  Control 
Unit  (Attenuator),  manufactured  by  the  David  Taylor  Model  Basin0   No  ampli- 
fication was  possible  with  this-- unit,  the  maximum  output  being  one  times 
the  input.  This  unit  also  permitted  a  convenient  switching  arrangement 
for  obtaining  a  calibration  trace  on  the  oscillograph. 

The  output  from  the  Control  Unit  was  picked  up  by  a  Consolidated 
Electrodynamics  Corporation  Recording  Oscillograph.   The  oscillograph 
utilized  galvanometers  of  the  fluid  damped  type  for  recording  the  velocity 
response  on  the  oscillograph  record.   The  galvanometers  had  a  natural 
frequency  of  1000  cps. 

Outputs  from  the  accelerometers  were  amplified  by  Consolidated 
Engineering  Corporation  3-KC  Carrier  Amplifiers.   The  power  supply  for 
the  above  amplifiers  consisted  of  a  Consolidated  Engineering  Corporation 
Oscillator-Power  Supply.   The  amplified  accelerometer  signals  were  then 
received  by  the  above  mentioned  oscillograph. 

A  General  Radio  Company,  1000  cycle,  Vacuum-Tube  Fork  was  used  to 
produce  a  time  reference  trace  on  the  oscillograph  record,,   The  oscillo- 
graph also  produced  a  time  reference;  however,  it  was  felt  that  the  1000 


cycle  trace  would  provide  a  good  checl  illograph's  time  sea 

Since  difficulty  was  encountered  with  the  interference  between  the 
test  structure  and  the  anvil  supporting  wires,  an  extra  galvanometer  was 
utilized  for  recording  an  indication  of  contact  between  the  wires  and  the 
structure.   This  consisted  of  a  simple  battery  circuit,  using  the  point  of 
probable  contact  as  the  switch.   When  wires  and  structure  made  contact,  a 
pulse  was  recorded  on  the  oscillograph  record,  the  duration  of  this  pulse 
gave  the  period  of  contact.   Figure  (4)  shows  a  sketch  of  the  instrumenta- 
tion. 

Appendix  VI  gives  a  tabulation  of  the  various  instrument  serial  and 
model  numbers. 


Test  Procedure 

The  following  steps  were  followed  for  a  typical  test  rum 

1.  The  structure's  columns  (springs)  were  placed  in  position  and  securely 
fastened, 

2.  The  anvil  was  elevated  to  a  position  where  interference  between  the 
structure  and  supporting  wires  would  not  be  encountered. 

3.  The  velocity  meter  magnets  were  adjusted  to  insure  free  passage  as  the 
anvil  moved  away. 

4.  The  hammer  was  elevated  to  a  predetermined  position,  depending  upon 
the  step  input  desired,  using  the  electric  hoist. 

5.  The  recording  oscillograph  was  operated  at  a  paper  speed  of  4  in. /sec. 
and  velocity  meter  and  accelerdmeter  calibration  traces  were  recorded. 

6.  All  calibration  settings  for  the  various  instruments  were  recorded. 

7.  The  instruments  were  then  set  for  the  predetermined  settings  for  the 
actual  run,  and  all  values  recorded. 

8.  The  recording  oscillograph  was  operated  at  a  paper  speed  of  40  in. /sec. 
and  the  hammer  released. 

9.  The  oscillograph  was  secured  when  the  velocity  meter  magnets  were  free 
from  the  coils  or,  for  the  lower  input  runs,  when  it  was  felt  sufficient 
data  were  recorded. 


List  of  Symbols 

a  -  initial  horizontal  anvil  displacement 

b  -  width  of  rectangular  supporting  column 

E  -  modulus  of  elasticity 

F.  -  force  produced  by  flexure  of  columns  (i  -  1,2) 

2 

g  -  acceleration  of  gravity  (32.2  ft/sec  ) 

h  -  thickness  of  rectangular  supporting  columns 

H  -  (subscript)  denotes  horizontal  component 

I  -  moment  of  inertia  of  the  cross  section 

k  -  spring  constant  in  elastic  range 

k*  -  spring  constant  in  plastic  range 

J^  -  column  length 

L  -  length  of  anvil  supporting  cable 

M.  -  mass  (i  =  1,2,3) 

M  -  fully  plastic  bending  moment 

Q  -  axial  load  on  columns  produced  by  structures  weight 

R  -  tension  in  hold-back  line 

T  -  tension  in  supporting  cables 

V  -  (subscript)  denotes  vertical  component 

x.  -  horizontal  displacement  of  M. 
l  i 

x  -  relative  spring  displacement  at  which  the  fully  plastic  bending 
moment  is  realized 

y  -  vertical  displacement  of  anvil 

(T  -  dynamic  yield  strength 
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3.   Theoretical  Considerations 
Mathematical  Model 

The  mathematical  model  of  the  two-story  structure  can  be  visualized 
as  a  mass  -  spring  system  as  shown  in  figure  (5a).   This  idealization 
reduces  the  structure  to  a  three  degree  of  freedom  system.   Figure  (5b) 
shows  the  free  body  diagram  for  the  system.   The  M.  symbols  denote  the 
masses  of  the  anvil,  center  mass,  and  upper  mass  respectively,  the  F. 
symbols  denote  the  forces  produced  by  the  flexure  of  the  columns,  and 
the  X.  denote  displacements  of  the  respective  masses. 

Summing  forces  for  each  mass  results  in  the  following  equations  of 
motion, 

(1) 

(2) 

-«  -       „5  <3) 

where  the  force  T   is  the  horizontal  component  of  the  force  in  the 
supporting  cables. 

Letting  ML"  denote  the  length  of  the  anvil  supporting  wire,  "a" 
the  horizontal  distance  that  the  anvil  is  initially  displaced,  "x',',   the 
displacement  of  the  anvil  from  the  initial  position  "a",  and  "vi",  the 
vertical  displacement  of  the  anvil,  we  have  from  figure  (6a), 

Lf  =  (*-  +  x,)a  +  C^-Sif 

where,  ^  =  YL2-  a.* 

In  figure  (6b)  "W"  is  the  total  weight  of  the  anvil  plus  structure 
(W  =  656  lbs.).   The  symbol  "T"  denotes  the  tension  in  the  supporting  cables, 
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*,    =  - 

M, 

•  * 

x2    = 

F,-Fz 
M2 

•  • 

y      — 

R 

The  symbol  "R"  denotes  the  hold-back  wire  reaction  which  immediately 
reduces  to  zero  upon  impact.   From  figure  (6b)  we  may  write, 

Tv  =   65G  +   -&!&■  ^  (4) 


Tv      Coe>  -Q-  ^a y. 


(5) 


from  figure  (6a), 


¥,  =  £    -  Vl2-  (a-f-X,)2  (6) 

V  -     C^+M  X,  m 
'  ~~      £  -  2/,  l  J 

V  —    (o.->-x,x)  X,  -r-  X,  +-  y.  /gx 

By  substitution  of  equation  (8)  into  equation  (4)  and  then  combining 

with  equation  (5),  we  arrive  at  a  value  for  T  in  terms  of  known  or 

H 

computed  variables.   Substituting  this  value  of  T  back  into  equation 

H 


(1)  and  solving  for  X   produces, 


3 


The  three  equations  describing  the  response  of  the  model  are  equa- 
tions (2),  (3),  and  (9). 

Basic  Considerations 

The  test  structure  followed  closely  to  the  portal  frame  used  by 
Kurzenhauser  in  his  tests,  the  major  differences  being  in  the  number  of 
degrees  of  freedom  and  the  method  of  loading.   The  portal  frame  had  a 
fixed  base  and  the  loading  took  place  through  the  upper  mass  of  the  struc- 
ture.  The  two-story  structure  used  in  these  tests  was  loaded  through  its 
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base. 

when  beycr       '  rate  red,  it  be- 

came necessary  to  deviate  from  Hooke's  Law  and  to  study  the  behavior  in 
the  plastic  region.   Since  the  entire  system  was  reduced  to  a  spring  - 
mass  system  it  became  necessary  to  assume  a  force  versus  deflection  relation- 
ship that  would  satisfy  the  behavior  of  the  structure  as  it  passed  from  the 
elastic  region  into  the  plastic  region.   Figure  (7)  shows  the  assumed 
force  versus  relative  deflection  response. 

The  structure's  columns  behaved  elastically  up  to  the  relative  deflec- 
tion corresponding  to  X*  ,  whereupon  the  behavior  departed  from  Hooke's  Law. 
The  simplest  way  to  represent  the  complicated  response  in  the  elastic-plastic 
deformation,  is  tc        two  linear  spring  constants,  (slope  of  the  lines 
represented  by  the  force  deflection  curve)  one  for  the  elastic  range  and  a 
second  for  the  plsstic  range.   The  elastic  spring  constant  is  represented 
by  k  and  the  plastic  spring  constant  by  k*.   Upon  loading,  the  columns 
follow  the  elastic  spring  constant  (k)  until  yielding  of  the  columns  takes 
place,  then  they  follow  the  plastic  spring  constant  (k*).  When  unloaded, 
the  columns  follow  the  displaced  elastic  spring  constant  (k),  forming  an 
hysteresis  loop  as  indicated  in  the  figure. 

It  is  assumed  that  the  transition  from  the  elastic  bendinj  to 
the  fully  plastic  bending  moment  in  the  column  cross-sections,  occurs  in- 
stantaneously at  the  deflection  corresponding  to  Xu., 

For  rectangular  cross-sections,  the  fully  plastic  bending  moment  is 

defined  by  the  relation,  z 

—  crb  h 

where  cP  is  the  yield  point  or  yield  strength  of  the  material  involved, 
b  is  the  width  of  the  column,  and  h  is  the  thickness  of  the  column. 
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A  force  analysis  of  one  of  the  columns  in  the  elastic  region  will 
permit  evaluation  of  the  elastic  spring  constant.   Referring  to  figure 
(8a),  and  utilizing  the  cantilever  beam  deflection  equations  gives  the 
following  results, 

r* 


jl»  ®  (4-) 

2       3  E  I  Where 


24-EI 


F  =  total  forcing  function 
on  twn  supporting 
colunns 

E  *  modulus  of  elasticity 
I  =  moment  of  inertia 


Considering  the  influence  of  the  compressive  force  Q,  where  Q  denotes  the 
total  compressive  force  on  two  columns  and  Q  is  the  critical  Euler  column 
loading  for  two  columns,  we  have, 


whe  re , 


K  = 


Q, 


_F 
X 


24E1 


ZU^EI 


(10) 


It  has  been  shown  that  the  approximation  used  in  equation  (10, 
sufficiently  accurate  almost  up  to  the  critical  value  of  the  load.  "   [6_| 

A  force  analysis  of  one  of  the  supporting  columns,  figure  (8b),  when 
M  is  achieved,  gives  the  following  for  small  deflections, 

SMa  =    o 
F  =-fx+  x  (H) 


Equation  (11)  is  the  equation  of  the  straight  line  representing  the 
plastic  region  shown  on  the  force-deflection  diagram  of  figure  (7).   The 
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spring  constant  k*  for  the  plastic  region  is  therefore  the  slope  of  this 
line, 

*  Q. 

K  =  -  -j-  d2) 

It  will  be  noted  from  figure  (7}  that,  at  the  deflection  x  ,  the 
spring  force  represented  by  equation  (10)  for  the  elastic  region  is  equal 
to  the  force  represented  by  equation  (12)  for  the  plastic  region.   Equat- 
ing these  two  produces  the  following, 


24  EI  Xu,  T|  _  Q. 


-  _JLX*  fA|k  (13) 


The  value  for  x  can  be  ea=sily  obtained  from  equation  (13). 


Complications  Introduced  by  Experiments!  Set-up 

The  two  major  problems  that  arose  from  the  experimental  model  were 
the  failure  to  record  the  initial  anvil  displacement  prior  to  impact  and 
the  lack  of  fixity  of  the  supporting  columns. 

Since  the  initial  anvil  displacement  was  not  recorded,  it  became 
necessary  to  determine  this  value  from  the  available  information   This  was 
accomplished  by  plotting  the  sum  of  the  system's  total  momentum  versus  time 
for  both  the  experimental  and  theoretical  values.  The  initial  displace-  , 
ment  "a"  used  in  equation  (9)  was  varied  through  a  reasonable  range  unti1 
agreement  with  the  experimental  points  was  obtained.   Figure  (9)  is  a 
representative  plot.   The  above  mentioned  comparison  between  theoretical 
and  experimental  was  made  for  each  run  where  sufficient  values  of  anvil 
velocity  were  available.   Since  the  length  of  the  velocity  trace  depended 
upon  the  time  the  magnet  departed  the  velocity  meter  coil,  tbe  values  for 
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the  anvil  velocity  for  the  plastic  runs  were  not  available  beyond  about 
50  to  1GC     (millis  onds).   The  value  obtained  for  the  initial  anvil 
displacement  utilizing  the  above  method  was  consistently  0.7  ft.   The 
value  arrived  at  for  this  displacement  should  have  been  fixed  for  all 
runs,  as  indicated  from  the  plots,  since  the  anvil  was  moved  to  the  same 
position  for  each  run  during  the  actual  tests. 

The  value  (a  =  0.7')  was  therefore  used  for  all  runs  with  the  ex- 
ception of  two  of  the  low  input  steel  tuns,  which  definitely  indicated  a 
lower  value  for  "a". 

The  lack  of  fixity  of  the  supporting  columns  introduced  a  problem  as 
to  the  value  for  effective  length  in  the  formulae  for  determing  k  and  k*. 
The  effective  length  was  varied  until  phase  agreement  (peaking  at  the  same 
time)  was  achieved.   The  effective  lengths  found  necessary  were  9,0  inches 
for  steel  and  8.75  inches  for  aluminum. 
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4.   Numerical  Solution 

Reduction  of  Equations  for  Computer 

The  first  step  in  reducing  the  equations  to  a  suitable  form  for 
introduction  into  the  computer  was  to  reduce  the  three  second  order 
equations  (2),  (3),  and  (9)  to  six  first  order  equations.   Introducing 
the  following  dummy  variables, 

•f  -  X, 

ft  =   x, 

y  =  x3 

then, 


J 

—~ 

X, 

n 

— 

•     • 

• 

zz 

•     • 

x3 

the  six  first  order  equations  then  become, 

X,  =    1  <14> 

X2  =  V  (15) 


x3=  *  <16> 


y  = 


(17)* 


3 


F,   -  Fa 


?  -     r'  ~  ^a  (18) 


*  =     -*-  (») 


*  Note  that  3^—  jli,  and  ^,  are  given  by  algebraic  relations  involving 
X|   and  $    i    see  Equations  (6)  and  (7),  page   12 
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Letting  the  symbol  D  denote  the  dependent  variable,  Y  denote  the 
independent  variable,  and  Introducing  the  following  equivalent  symbols, 

D,  =    S  Y,    =  X,  CC  =   a-t- X, 

IX  =    V  Y2  =  X,,  CK  =  £.-  ^, 

D,  =    y  Y3  =  X3  YV  =    Sf, 

D«.  =    J  Y^  =  :f 

D,  =7  Yff  =  ? 

d6  »   *  Y4  =  Y 

produces  the  following  equations  written  in  the  Fortran  language, 

DC)  =  Y(+) 
DW  =  YU) 
DCi)  =   Y(6) 

0(4)  =  (-F0)  -  Ccc/c/c)  *  656.  *  (/.  +  (  Y(*)**2  +  YV**-2)/ 
(  52.20  *00))/f/8. 3+656.  /32.20)*CCC/C/<)  **a) 

D(5)=   (F(0-  F(2))/o.7// 

DU>  -     F(2)  /0.705T 

Appendix  I  contains  the  Main  Flow  Chart,  the  Force  versus  Relative  De- 
flection Flow  Chart,  and  the  actual  Fortran  program.   The  program  traces 
the  time  history  of  the  displacements,  velocities,  and  accelerations* 
of  the  three  masses  M, ,  M0,  and  M  starting  with  the  initial  conditions 
X][  =  x2  =  x.,  =  0;  V2  =  V3  =  0,  Vx  =  STEP. 


^Accelerations  are  given  in  the  actual  computer  print-out,  he  river,  they  were 
not  used  in  the  analysis  and  have  been  trimmed  off  the  sheets  shown  in  Appen- 
dix V.   Also,  the  print  statement  shown  in  Appendix  I  does  not  produce  dis- 
placements.  Displacements  were,  however,  printed  in  an  earlier  version  and 
it  would  be  very  easy  to  restore  them. 
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Numerical  Integration  System 

The  numerical  integration  system  used  w        mgi  nod. 

The  two  major  advantages  of  this  system,  as  stated  by  Milne  |_PJ  ,  are: 

1.  No  special  starting  procedure  is  required. 

2.  The  length  of  the  time  interval  may  be  altered  during  the 
computation. 

The  reasons  for  utilizing  the  Runge-Kutta  method  are;:   fiiet  of  all  the 

fact  that  the  integration  relies  upon  the  last  computed  point,  as  listed 

in  advantage  number  one  above,  secondly,  the  routine  was  already  available 

for  immediate  use  in  the  U.  S.  Naval  Postgraduate  School  computer  center. 

The  system,  as  with  many  numerical  integration  systems,  is  based  on 

the  Taylor  Series.   For  a  single  first  order  differential  equation  y  - 

f (  "t  »y)   the  following  computations  are  completed  to  arrive  at  a  value 

y(  t  +  At), 

K,     =      Atf  (t^) 

K5     =•    Atfft^fj  ^J  +  i  Kz) 
K4    =     A-tf  U-h^-j  }-t-   Kz) 

The  same  basic  procedure  as  above  can  be  applied  to  a  system  of  first 
order  equations  such  as  the  six  equations  shown  previously. 
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This  section  shows  a  graphical  comparison  of  the  experimental  results, 
represented  by  the  cin  I   points,  and  the  theoretical  response  predicted 
by  the  mathematical  model,  represented  by  the  smooth  curves. 

The  first  section  displays  the  results  of  sixteen  runs  using  steel 
columns.   The  first  seven  plots  contain  a  single  theoretical  curve  since 
these  columns  experienced  no  plastic  deformation.   The  next  nine  plots  dis- 
play three  theoretical  curves  corresponding  to  dynamic  yield  strengths  of 
40,000  psi,  50,000  psi,  and  60,000  psi.   The  last  two  plots  repeat  two  of 
the  runs  extending  the  time  cut-off  to  200  MS. 

The  second  section  displays  the  results  of  ten  runs  using  aluminum 
columns.   The  first  three  runs  arja  in  the  elastic  region  and  therefore 
display  the  single  theoretical  curve.   The  remaining  seven  runs  display 
two  theoretical  curves  for  dynamic  yield  strengths  of  20,000  psi  and  30,000 
psi.   The  theoretical  curves  shown  on  pages  44  ,  45  ,  and  46  actually 
merge,  within  plotting  error,  from  approximately  90  MS  to  100  MS. 

The  term  "Step"  denotes  the  initial  anvil  velocity  upon  impact  (ft/ 
sec).   It  was  determined  from  the  velocity  meter  mounted  on  the  anvil. 
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Section  One  -  Steel 
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Section  Two  -  Aluminum  Alloy 
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6.   Discussion  of  Results 

1.  A  complete  tabulation  of  the  eyperimental  and  theoretical  points  shown 
in  the  graphical  comparison  is  not  included  in  this  thesis.   Samples  of 

the  data  are  contained  in  Appendices  (II)  and  (V).   The  Fortran  program  used 
is  found  in  Appendix  (I),  and,  if  desired,  may  be  utilized  for  obtaining 
additional  theoretical  results.   The  experimental  results  shown  on  the 
various  plots  are  as  accurate  as  the  values  used  in  plotting. 

2.  The  two  major  experimental  discrepancies,  as  noted  earlier,  were  the 
lack  of  fixity  and  the  failure  to  record  the  initial  anvil  displacement. 
These  two  items  have  been  suitably  discussed  in  previous  sections.;  however, 
an  improvement  to  the  test  structure  is  suggested  for  future  tests.   The 
single  bolting  arrangement  used  appeared  to  be  inadequate,  as  far  as  pro- 
viding a  fixity  approaching  a  cantilever  beam.   The  more  rigid  connection 
used  by  Kurzenhauser  [_1J  in  his  tests  consisted  of  four  5/16"  N.  C.  socket 
head  cap  screws  plus  a  1/4"  clamping  plate.   The  four  screws  were  located 
near  the  edges  of  the  clamping  plate  instead  of  the  centrally  located  bolt 
used  in  this  thesis.   It  is  therefore  felt  that  the  design  of  the  test 
structure  was  poor  in  this  respect,  and  the  use  of  four  holding  bolts 
would  have  improved  the  fixity. 

Nothing  further  will  be  said  concerning  the  failure  to  record  the  in- 
itial anvil  displacement,  except  that  this  was  vital  information  that  should 
have  been  recorded.   In  any  future  structural  tests  utilizing  the  ballistic 
pendulum  herein  described,  the  test  structure  should  receive  the  initial  im- 
pact with  the  anvil  hanging  vertically  (initial  displacement  zero),  thus 
eliminating  this  initial  horizontal  force. 
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3.   The  maximum  step  velocity  (Initial  anvil  velocity  upon  impact)  ex- 
perienced by  the  test  structure  was  7.83  ft/sec.   The  ballistic  pendulum 
was  capable  of  producing  approximately  14.0  ft/sec;  however,  due  to  the 
interference  between  the  test  structure  and  the  anvil  supporting  cables, 
the  maximum  of  7.83  ft/sec.  was  the  highest  input  before  interference  was 
experienced. 

4„   The  limiting  value  of  time  used  for  the  graphical  comparison  of  results 
was  100  MS.   The  reason  for  cutting  off  the  comparison  at  this  point  was 
due  to  the  fact  that  in  reconstructing  the  initial  anvil  displacement 
(a  =  0.7'),  the  total  momentum  vs.  time  plot  was  used.   This  plot  was 
limited  by  the  velocity  record  of  the  anvil.  The  velocity  meter  magnet 
for  the  anvil  departed  the  coil  after  approximately  100  MS,  up  to  a  step  of 
3.86  ft/sec.   Above  this  value  for  initial  velocity,  the  anvil  velocity 
trace  was  even  shorter,  averaging  about  50  MS.   Since  the  total  momentum 
vs.  time  plot  was  limited  to  approximately  100  MS,  it  was  felt  that  use  of 
the  predicted  anvil  displacement  beyond  100  MS  would  not  be  warrented. 

A  second  reason  for  using  the  100  MS  limit  was  that  the  two  velocity 
meters  mounted  on  the  structural  masses  produced  faulty  traces  beyond  70  MS 
for  the  steel  columns,  at  the  higher  inputs.  These  faulty  traces  were  like- 
wise caused  by  the  velocity  meter  magnets  departing  their  respective  coils. 
The  purpose  of  mounting  the  accelerometers  on  the  test  structure  was  to 
allow  for  a  continuation  of  the  velocity  plots  beyond  70  MS  by  integration  of 
the  acceleration  vs.  time  traces.   It  was  felt  that  the  increased  error  due 
to  reduction  and  integration  of  the  acceleration  traces  would  be  held  to  a 
minimum  if  the  number  of  points  added  by  this  method  were  small.   Three 
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points  were  added,  producing  a  total  velocity  vs.  time  trace  of  100  MS 
duration.   Discrepancies  of  the  order  i  5%  were  noted  between  velocities 
obtained  from  the  velocity  traces  directly  and  by  integration  of  the 
acceleration  traces,  over  intervals  where  the  velocity  traces  were  consider- 
ed valid. 

5.  The  graphical  comparison  of  results,  for  steel  columns,  found  on  pages 
22  through  37  ,  consists  of  approximately   seven  elastic  runs  and  nine 
plastic  runs.   The  experimental  points  for  the  elastic  runs  are  compared 
with  a  single  theoretical  plot  since  the  theoretical  tesults  within  this 
range,  are  not  affected  by  changes  in  the  dynamic  yield  strength.   The  plot 
for  step  =  2.87  ft/sec  is  a  borderline  case.   A  slight  variation  between 
theoretical  values  for  cr  =  40,000  psi  and  c3  =  60,000  psi  was  noted,  how- 
ever the  variation  was  quite  small  (0.3  ft/sec  maximum)  and  therefore  only 
the  theoretical  trace  f or  ^ =  40,000  psi  is  shown.   Starting  with  the  value 
step  «  3.86  ft/sec,  three  theoretical  curves  are  shown  for  the  remainder 
of  the  runs.   The  curves  represent  dynamic  yield  strengths  of  40,000  psi, 
50,000  psi,  and  60,000  psi.   The  general  trend  in  the  curves  indicates  the 
increased  dynamic  yield  strength  with  increased  initial  velocity.   The  runs 
for  step  =  5.02  ft/sec  and  7.23  ft/sec  do  not  fit  in  with  the  remainder  of 
the  plots.   The  dynamic  yield  strength  for  these  two  runs  shows  a  decrease, 
whereas,  an  increase  should  be  indicated  to  follow  the  sequence  of  values 
previous  to  and  following  these  respective  runs.   The  cause  of  the  dis- 
agreement could  possibly  be  due  to  an  error  in  the  experimental  values  for 
initial  anvil  velocity.   The  comparison  indicates  a  lower  velocity  input 
than  the  values  used.   A  change  of  0.02"  in  measurement  from  the  oscillo- 
graph record  produces  a  change  of  0.19  ft/sec.   Also  a  variation  in  support- 
ing columns  from  one  test  structure  to  another  could  introduce  an  appreciable 


error.   It  can  be  seen  in  some  of  the  higher  input  plots,  for  steel,  that 
the  experimental  points  tend  to  fall  off      rapidly  than  would  be  indicat- 
ed for  continuity  of  plot.   As  w   acted  earlier,  the  last  three  points 
for  the  steel  plastic  runs  were  computed  using  the  integrated  acceleration 
vs.  time  t races  .   The  added  error  caused  by  this  procedure  may  account  for 
this  discrepancy. 

Pages  38  and  39  show  an  extension  of  two  of  the  previous  runs  out  to 
200  MS.   The  purpose  of  these  two  plots  was  to  give  an  indication  of  the 
agreement  beyond  the  100  MS  limit.   As  can  be  seen  from  the  plots,  the 
agreement  is  very  good.   This  further  indicates  the  validity  of  the  computed 
value  for  initial  anvil  displacement. 

- 

6.  The  graphical  comparisons  using  aluminum  columns  are  likewise  limited 

to  100  MS.   The  first  three  plots  show  the  single  theoretical  trace  ivhereas 
the  remaining  seven  plots  show  two  theoretical  traces,  one  trace  representing 
CT     ?0,000  psi  and  a  second  trace  representing  CT  =  30,000  psi.   For 
the  three  runs  represented  by  step  -  2.29,  2.49,  and  2.82  £t/sec,  the 
theoretical  values  for  ifr    =  20,000  psi  and  cr  .  30,000  psi  merge  within 
plotting  error  between  90  and  100  MS. 

The  experimental  values  for  the  velocity-time  plot  for  the  aluminum 
columns  were  taken  entirely  from  the  velocity  trace  since  the  magnets  re- 
mained in  the  coils  somewhat  longer  than  for  the  steel  columns.   For  the 
two  highest  input  velocities  the  last  three  points  were  checked  by  integrat- 
ing the  acceleration  vs.  time  traces. 

7.  Massard  and  Collins  [2J  ,  in  tests  of  structural  metals  under  rapid 
loading,  demonstrated  the  time  sensitive  behavior  known  as  "delayed  yield" 
for  structural  steel  and  aluminum.   It  was  found  that  metals  that  display 
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a  definite  upper  and  lower  yield  point,  under  slow  testing  procedures, 
are  quite  time  sensitive  and  display  the  delayed  yield  phenomenon.   For 
metals  such  as  aluminum,  where  straining  under  slow  tests  is  not  dis- 
continuous, the  delayed  yield  effect  is  not  as  prominent.   From  the 
graphical  presentation  of  data,  the  aluminum  plots  would  tend  to  substan- 
tiate this  conclusion,  since  the  apparent  increase  in  dynamic  yield 
strength  is  somewhat  smaller  relative  to  that  demonstrated  for  steel. 

8.  An  improvement  to  the  testing  procedure  used  would  be  to  replace  the 
supporting  columns  before  each  run.   Obviously,  this  is  a  requirement  for 
the  plastic  runs;  however,  in  several  instances,  the  same  set  of  supporting 
columns  was  used  for  one  very  low  input  elastic  run  and  then  a  higher 
velocity  input  plastic  run.   This  procedure  only  involved  about  three  of 
the  runs  and  was  done  to  save  time. 

9.  As  mentioned  earlier,  the  integration  routine  utilized  in  this  thesis 
was  the  Runge-Kutta  method.   The  time  interval  used  was  1/10  MS.   This  was 
arrived  at  by  comparison  with  the  results  for  an  exact  solution  for  a 
similar  system  contained  in  reference  [j4]  .   The  only  difference  between 
the  two  systems  was  that  the  exact  solution  was  based  upon  a  system  whose 
base  was  rigid.   The  system  utilized  in  this  thesis  was  easily  reduced  to 
the  above  system  merely  by  making  M.  very  large.   Comparison  of  the  two 
systems  in  the  elastic  range  (exact  solution  was  valid  for  elastic  response 
only)  produced  answers  that  were  in  agreement  to  about  three  significant 
figures.   This  was  considered  accurate  enough  for  plotting. 

Since  the  time  interval  used  in  the  computer  program  was  1/10  MS, 
theoretical  results  could  have  been  obtained  for  every  1/10  MS;  however, 
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this  would  have  produced  an  excessive  amount  of  output  data  and  the  values 
were  therefore  printed  every  MS  or  every  tenth  time  interval.   The  theo- 
retical curves  shown  in  the  graphical  presentation  of  results  were  plotted 
using  points  every  five  milliseconds. 

One  inherent  difficulty  with  the  Runge-Kutta  method,  as  pointed  out 
by  Newmark  and  Chan  [si  ,  is  that  "the  amplitude  will  gradually  damp  it- 
self out  even  in  the  case  of  an  undamped  system,  and  the  method  is  there- 
fore not  desirable  for  a  long  period  of  time."  No  explanation  was  given  as 
to  what  constituted  "a  long  period  of  time";  however,  the  time  limit  imposed 
in  this  problem  was  undoubtedly  short  enough  to  prevent  this  from  affecting 
the  results. 

10.   Perhaps  the  greatest  source  of  error  involved  in  the  experiment  was 
the  reduction  of  the  experimental  data.   The  data  reduction  involved  first 
of  all  a  determination  of  a  zero  reference  level  for  the  various  traces. 
This  was  accomplished  by  using  the  edge  of  the  oscillograph  record  as  a 
reference  and  then  comparing  this  with  the  1000  cycle  trace  for  straight- 
ness.   In  all  cases  the  straight  line  formed  by  the  1000  cycle  trace  agreed 
with  the  paper  edge  reference.   Amplitudes  for  the  various  parameters  were 
measured  from  the  respective  reference  lines  using  a  six  inch  steel  scale. 
A  magnifying  glass  was  used  for  greater  accuracy.   The  scale  was  read  to 
the  nearest  0.01". 

As  noted  earlier,  the  sensitivity  of  the  velocity  meters  varied  ±  37=, 
during  a  4"  sweep.   The  response  of  the  accelerometers  was  good;  the  only 
drawback  in  their  use  was  the  high  frequency  pick-up  in  the  first  30  MS. 
This  completely  obliterated  the  trace  for  the  first  30  or  40  MS,  prevent- 
ing any  real  comparison  of  the  compatibility  between  velocity  and  accelera- 
tion traces. 
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The  intec  the  aecelc     i  vs.  time  plots  was  accomplished 

by  using  a  K  &  E  4236  Compensating  Polar  Planimeter,  serial  number  6098. 
It  is  estimated  that  an  error  of  about  2%  was  possible  using  this  instru- 
ment since  the  areas  involved  were  of  the  order  0.1  sq.  in. 

Since  the  values  obtained  in  the  data  reduction  were  accurate  to  about 
two  or  three  significant  figures,  the  values  put  into  the  computer  were 
limited  to  the  same  precision,  which  was  also  utilized  in  presenting  the 
graphical  comparison  of  results. 

11.   Tt  is  to  be  noted  that  the  parameters  chosen  for  the  comparison  of  the 
results  could  have  been  force  vs.  deflection,  acceleration  vs.  time,  or  some 
other  comparison.  The  velocity  ,vs.  time  comparison  was  chosen  because  it 
offered  the  quickest  comparison  with  the  least  manipulation  of  numbers.  The 
other  comparisons  mentioned  above,  would  have  involved  extensive  reduction, 
which  would  have  introduced  the  possibility  of  increased  error. 
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7.   Conclusions  and  Recommendations 

The  gratifyingly  close  comparisons  between  e xperimentally  observed 
behavior  and  theoretically  predicted  behavior  of  the  three  mass  structure 
considered  herein  permits  the  conclusion  that  the  response  of  si 
structures  may  be  predicted  by  mathematical  analysis  using  digital  computer 
calculations  and  involving  force-deflection  relations  for  the.  elastic  and 
plastic  regimes  which  Im  '    prior  deformation  history  in  the  manner  des- 
cribed in  Appendix  T.   The  analysis  in  the  prese      c  was  complicated  by 
a  feature  of  the  experimental  set-up  which  would  probably  have  no  counter- 
part in  actual  practice.  Tl      oald  be  no  ess     '  difficulty  in  ana ! 
ing  cases  where  the  excitation  comes  from  lateral  forces  havinj      i  f  ied 
variations  with  time  rather  th         velocity  step  of  one  ses, 

although  the  latter  case  is        ly  of  primary  interest  in  considering 
response  to  explosive  loadings   Similarly,  there  should  be  no      culty 
in  incorporating  further  refinements  regardii    I  i ge"  geometi      he 
analysis  should  it  be  desired  to   1T    the  deflection  history  int     .^es 
where  the  present  restrictions  to  "small"  geometry  might  n  ciently 

accurate. 

The  present  experi  t  for  the  -ities 

that  accompany  the  larger  input       ty  steps,  the       of  cr        es 
the  best  theoretical  prediction  of  the  respon;    I  o  Increases,,   For  rea; 
cited  by  Kurzenhauser  [ij  ,  no  attempt  was  made  in  this  thesis  : 
quantitative  evaluation  of  this  i      >n. 

The  predicted  response      Is  greatly  on  the  details  of  the  load 
deflection  relations  which  depend,  in  turn,  upon  a  knowledge  of  the 
behavior  of  the  material  as  it  is  loaded  and  unloaded  in  both  the  elastic 
and  plastic  regime.   It  woul        great  advv         have  better  informa- 
tion than  seems  to  be  available  j      ly  about  structural  materials  loa 
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in  such  a  manner;  including  the  variability  of  cr  with  strain-rate. 

Future  experimental  investigations  should  employ  apparatus  in  which 
a  closer  control  over  or  knowledge  of  column  fixity  is  available.   Also, 
if  a  ballistic  pendulum  type  of  apparatus  is  to  be  used,  it  is  clear  that 
the  masses  should  be  hanging  freely  before  impact  so  as  to  avoid  the  compli- 
cations introduced  in  the  present  experiments  by  the  presence  of  the  hold- 
back cable.   Finally,  it  would  clearly  be  quite  advantageous  to  have  a 
single  instrumentation  capable  of  reporting  the  response  of  the  structure 
beyond  the  100  MS  cut-off  without  the  necessity  of  augmenting  the  response 
with  a  second  instrumentation  as  was  necessary  in  this  investigation,  where' 
in  velocity  meter  data  was  extended  by  use  of  integrated  accelerometer 
data. 
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Figure  9 
Plot  of  total  momentum  vs.  time  for  an  initial  velocity  input 
(Step)  of  1.82  ft/sec.   The  symbol  O  denotes  experimental 
points  and  the  solid  curve  represents  theoretical  values  for 
an  initial  anvil  displacement  of  0.7  ft. 
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; ] ENDIX  I 
FLOW  CHARTS  AND  COMPUTER  PROGRAM 
This  section  contains  a  very  general  flow  chart  for  the  main 
program,  a  more  detailed  flow  chart  of  the  force  versus  deflection 
routine  which  was  the  controlling  routine  of  the  program,  and  the  actual 
computer  program.   The  force  versus  relative  deflection  routine  permitted 
the  evaluation  of  the  spring  forces  as  the  mathematical  model  experienced 
elastic  -  plastic  response.   A  list  of  symbols  is  included,  since  several 
symbol  changes  were  necessary  in  conversion  to  the  Fortran  language.   The 
figure  shown  belo\  and  che  force- deflection  flow  chart  sbown  on  page  74 
explain  the  force-deflection  subroutine. 


LOCUS     OF   MEW    CEM"i 


The  subroutine  is  er.ered  from  the  main  program  shown  by  the  box  labeled 
MP.   Several  parameters  are  picked  up  such  as  the  relative  displacement  of 
the  columns,  denoted  by  S,  the  abscissa  B  of  the  knee-point,  the  elastic, 
and  plastic  spring  constant  s  k  and  k*,  the  value  for  F  \;here  the  extension 
of  the  line  represented  by  k*  crosses  the  vertical  or  force  axis  (denoted 
by  A),  and  the  coordinates  (P,G)  of  the  present  location  of  the  origin  of 
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the  force-deflection  plot.   Note  that  k*  is  negative  as  derived  on  page  15 
This  is  indicated  in  the  figure  above.   The  symbol  R  is  the  value  of 
displacement  used  throughout  the  routine.   This  value  is  continually 
corrected  for  the  changing  origin.   The  location  of  the  origin  is  deter- 
mined from  the  previous  history  of  the  response.   The  first  time  that  the 
subroutine  is  entered  the  values  of  P  and  G  are  zero,  that  is  the  origin 
is  at  (0,0).   If  the  relative  displacement  S  is  less  than  the  knee-point 
value  B,  then  the  force  will  be  equal  to  the  elastic  spring  constant  k 
times  the  relative  displacement  R,  where  R  =  S  -  P  or  R  =  S  since  P  = 
zero.   When  the  value  of  R  builds  up  to  exceed  the  knee-point  value  B, 
a  check  is  made  to  see  if  the  column  is  loading  or  unloading.   The  check 
used  is  a  displacement  check  denoted  by  the  symbol  U,  where  U  is  equal 
to  the  present  relative  displacement  S  minus  the  previous  relative  dis- 
placement SS.   As  long  as  the  present  displacement  is  larger,  the  forcing 
function  is  equal  to  A  +  k*R.   The  above  routine  is  followed  until  the 
columns  begin  to  unload.   The  present  value  of  relative  displacement 
becomes  smaller  than  the  previous  value  and  a  new  origin  is  computed 
from  the  relations:      P  =  0  +  H  -  B  and  G  =  k*(P) 

where  0  is  the  horizontal  value  of  the  previous  origin  and  H  is  the  last 
recorded  value  of  R.   The  routine  loops  back  to  R  =  S  -  P  where  R  will 
now  be  less  than  the  knee-point  value  B.   The  value  for  force  will  be 
computed  from  F  =  G  +  kR  or  for  values  greater  than  BbyF=A+G+  k*R. 
The  above  described  routine  continues  until  the  problem  ends  and  a  new 
step  input  is  introduced  whereupon  all  parameters  are  initialized.   The 
same  type  of  routine  is  also  included  for  the  case  where  the  relative  dis- 
placement or  forcing  function  or  both  become  negative.   This  can  be  seen 
by  following  the  same  type  of  analysis  used  above. 
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The  force-deflection  routine  is  found  inside  the  subroutine  labeled 
SUBROUTINE  DERIV  statement  number  77  of  the  program  listing. 
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List  of  Symbols 

STEP  -   the  initial  input  velocity  to  anvil  (ft/sec) 

V    -  time  increment  (sec) 

XE   -  time  limit  of  problem  (sec) 

DE   -  end  value  of  step  input  (ft/sec) 

DD   -  initial  anvil  displacement  [a|*   (ft) 

CHG   -  correction  factor  to  allow  for  variable  initial  anvil 
displacement 


DC 


1    2    r  1     2   ._  2, 


i  anvil  cable  length]  '  -  [dd]     (ft  ) 


YD.   -  displacement  corresponding  to  ultimate  bending  moment 

2 
SL.   -  elastic  spring  constant  ^  [k]    (#/in  ) 


STR 


R.  -  plastic  spring  constant  [k*]   (#/in  ) 


FBR.  -  intersection  of  plastic  force-deflection  line  with  the 
vertical  axis  (lbs) 

ANS  -  output  printing  interval  (sec) 

SP(I)-  step  input  values  for  I  different  values  (ft/sec) 

Y.  -  independent  variables 

D.  -  dependent  variables 

F  -  forcing  function  in  columns  (lbs) 

S  -  present  relative  displacement 

SS  -  previous  relative  displacement 

X  -  time  variable 

TAB  -  time  expired  between  print  intervals 

0    -  last  value  for  relative  displacement  of  origin  of  force- 
deflection  plot 

H    -  last  value  for  force  origin  of  force-deflection  plot 


*Symbols  in  ["  I  refer  to  previously  defined  terms, 
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G    -  new  force  coordinate  for  force-deflect  ion  plot 

P    -  new  relative  displacement  coordinate  for  force-deflection 
plot 

R    -  Relative  deflection  of  the  columns 

U    -  S  — SS 

SUM  -  total  momentum  of  the  structure 
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Main  Flow  Chart 


Statement  // 


Yes 


11B1- 


Set  test  values, 
Time  interval  and 
Time  of  run 


Initialize 


Introduce  new  STEP 


Set  Anvil 
Velocity  =  STEP 


Augment  Time 


Evaluate  Force 
Functions  and 
Integrate  Eqns, 
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46 
53 
54 

55 

i 
62 
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**  Statement  numbers  are  sequential  starting  with 
PROGRAM  KENBOB  as  statement  number  one. 
Continued  statements  count  as  one  statement. 
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Computer  Program  for  Steel  Columns* 
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o— 

*The  program  for  aluminum  columns  is  identical  except  for 
slightly  different  values  for  M~  and  M_. 
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APPENDIX  II 
REDUCTION  OF  DATA 

This  section  contains  a  complete  set  of  calculations  for  one  run. 
Also  included  is  a  tracing  of  the  oscillograph  record  for  the  same  run. 
A  photograph  of  the  actual  record  was  not  possible  due  to  the  lack  of 
contrast  between  the  trace  and  the  paper.   It  is  to  be  noted  that  the 
1000  cycle  trace,  mentioned  earlier,  is  not  shown.   This  particular 
trace  was  omitted  to  prevent  cluttering  of  the  other  tracings  and  also 
because  of  the  tedious  tracing  job  required  to  accurately  reproduce  this 
trace. 

The  data  reduction  consisted  of  the  following  steps: 

1.  Establishing  reference  lines  on  oscillograph  record 

2.  Measuring  trace  amplitudes  for  velocity  parameters 

3.  Integrating  the  acceleration  vs.  time  trace 

4.  Computing  conversion  factors 

5.  Converting  measured  parameters  to  proper  units. 

The  run  in  question  was  for  steel  columns  and  a  step  input  of  5.02 
ft/sec. 

Computation  of  the  necessary  spring  constants,  and  other  input  para- 
meters are  also  shown. 
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Trace  of  Oscillograph  Record 
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Time 

V(l) 

V(2) 

0 

0.88 

0 

10 
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20 
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.59 

30 

* 

.75 

40 

* 

.75 

50 

* 

.86 

60 

* 

1.14 

70 

* 

1.39 

The  preceding  figure  shows  traces  of  the  oscillograph  record. 
The  symbols  V  denote  velocity  and  A  denote  acceleration.   Direct 
measurements  from  the  figure  gives  the  following  values  in  inches; 

V(3) 

0 

0 

.08 

.37 

.85 

1.24 

1.40 

1.37 

*values  not  measured 

Calibration  pulse  measurements: 

V(l)  -  0.720  inches 

V(2)  -  0.785  inches 

V(3)  -  0.820  inches 

A(2)  -  1.24  inches 

A(3)  -  1.12  inches 
Velocity  meter  sensitivities  and  accelerometer  calibration  constants 

V(l)  -  146  mv/fps 

V(2)  -  167  mv/fps 

V(3)  -  151  mv/fps 

**A(2)  -  70.7  g 

**A(3)  -  78.0  g 


**  values  obtained  by  the  David  Taylor  Model  Basin  using  calibration 
resistance. 
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Instrument  settings; 
a.   Calibration 


Amplifier  gain  set  at  30  for  accelerometers,  paper  speed 

°«t  at  4*'/sec 

Battery  voltage  for  velocity  meters  set  at  0.6  volts 

b.   Run 

Amplifier  gain  set  at  5  for  accelerometers,  paper  speed 
set  at  40"/sec 

Combination  of  the  above  information  gives  the  following  conversion 

factors: 

a.      for  velocity  conversion 

VfaCTOR     =  CAL1&.    VOLTAGE  Cv)         (  I 

S  EMS  I  Tl  V  IT   <  C  V/rps  )  PUL-S  £.      ht. 

SAMP; 


Vi  PACTOe     -       ^.5  7        //n. 


b.   for  acceleration  factors 


Af 


'AC  TDK      — 


-     CALie.  constant  run  cahm 


CAi  I  S       PUL5I  CALIB..  &  A  J    .0 


<\ 


*  mctor    -  1.24-     m.  To    "     ^°  /iec"' 

A5p^CTO(p     =       3  1  A-        /sec.-t»>. 

c.   conversion  for  integration  constant 

Vc-FAC-ra*  =  A  p«to«  *        TXME^CTDfc        •      AREA 

V2C .  FAc-yi>e       -         7.65"  •    /UJ//\      (  fV^-c  .) 
V3c>PrtcT»R  9.5  5    -AREA      C^V&eJ 
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Areas  under  acceleration  ttace 


Time  range  (MS) 
70=80 
70-90 
70-100 


Area  V(2)   (in  ) 

0.01 

-0.14 

-0.32 


Area  V(3)   (in  ) 
-0.04 
-0.06 
-0.14 


From  the  above  informations 


ime  (MS) 

VI 

:d 

(f 

0 

5.02 

10 

* 

20 

* 

30 

* 

40 

* 

50 

* 

60 

*  '■ 

70 

* 

80 

* 

90 

* 

100 

* 

V(2)   (ft/sec)    V(3)   (ft/sec) 


0 

0.82 

2.70 

3.43 

3.43 

3.92 

5.21 

6.35 

6.43 

5.28 

3.90 


*Not  measured 
Evaluation  of  spring  constants  -  steel 


0 

0 

0.39 

1.79 

4.12 

6.00 

6.78 

6.63 

6.26 

6.07 

5.32 


Column  thickness  =  h  =  0.128" 
Column  width     =  b  =  1.95" 
Effective  length  =  1    =  90." 

*modulus  of  elasticity  =  E  =  29.4  x  10  psi 

Weight  on  upper  columns  =  Q_  =  23.1  lbs 

Weight  on  lower  columns  =  Q  =47.0  lbs 


*Value  obtained  from  actual  tension  test  of  sample 
column.   See  Appendix  III 
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I  = 

0^- 


12  12 


2.1T'  2TTa'Ca9-4-)(5'gM) 


K,  = 


Jt3 


Q 


Ki 

K2 


b/i 


=     324-     ■  vin 
=i   52  7     ,b-/- 


K*  = 


41  0 


+  O  lbs 


i  1      =        (aMU9   4-)(54Q 

A  J  ^  3 


"l  -  ^2-1 


5.22       ,b/,». 


K. 


2.57      ,bS** 


Evaluation  of   force-deflection  knee  points 


(a..)  CJ"     —         4rQ  j   o  o  o      ps  «. 

2. 

•I-M^   =       ar  b  h/     -    C4-o,oooXi'3sX°  ,2©^    -     l27Si*-ib 


K ,  X  lA ,     — 


4- Ma   —    QiXhi 


32.4-  Xui   - 

X^  j  ~ 

Xo(2.  =■ 

A  - 


1218     —     47.0^, 


I2  7S 

19  6  3 

1273 

29  6  6 

<^   <Mw 

■sr        O.  4  3  I 


izlB 


je 


I  4 -2.       lb 
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( b )  cr   -  ,000    p s  c 
x^2  -_  0.53s      >  n 

A        -  '77       /I 

(c )  0s    =  GO,  000     ps  «. 
Xai    =■  O-  (6A-  6       /* 

A       =  2/3     lb. 

Using  the  same  procedure  for  3003-H14  aluminum: 

K,     =  I  03      '  b7ih  WH 


'/ 


b  -     I.9S 

K2     -       M2     ,b'/.*.  J, -.. 

*  ,-  fa 

t    =    /  o,  o  x   /  0   f>: 


lb 


K,       -     -    5.2  I     'Via,  Q(   :     22,  8     Ibj 


K2*    =    -2.53  'V,, 

KMtE      PoirJT£    1 

CO  <T"      -       20,000  ps  c 

Xu,    -      0.  fc  f  3  ;  v, 

X  u  t   -.:      0,616  '  ^ 

A        -      70,  8  lb. 

(b)        cr~     -     $o,ood  j 

Xu,     _       0.^28  IH, 

X «  0  .  <7  2  5  /  *  . 

A         -10  6  lb. 


Q-,--     /  /  .  ^    ■ 


*value   taken   from  Alcoa  Aluminum  Handbook   for   3003-H14  Aluminum 
reference      j  7|     . 
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APPENDIX  III 
COLUMN  MATERIALS,  TENSION  TESTS,  AND  EFFECTIVE  MASSES 
Column  Materials 

As  noted  earlier,  two  different  materials  were  utilized  in  the  test 
coluirns;  hot  rolled  steel,  with  properties  similar  to  those  of  the  steel 
used  by  Kurzenhauser,  and  3003  -  H14  aluminum.   The  aluminum  sheet  used 
had  a  yield  strength  of  21,000  psi  and  an  ultimate  strength  of  22,000 
psi,  as  listed  by  the  Alcoa  Aluminum  Handbook  [7J  .   A  6061  T6  aluminum 
used  by  Massard  and  Collins  |  2J   ,  would  have  been  more  representative 
of  the  aluminum  used  in  structural  work;  however,  this  was  not  available 
at  the  time  the  experimentation  was  conducted. 

Tension  Tests 

Three  samples  of  each  of  the  column  materials  used  were  selected  and 
shaped  into  standard  tensile  test  specimens  with  a  2"  gage  length  to  meet 
ASTM  Standard  E8-54T.  A  Baldwin  Universal  Testing  machine,  along  with  a 
Riehle  extensometer  were  utilized  in  making  the  tension  tests.  The  results 
are  tabulated  in  Table  1  below: 


Table  1 


Steel  Aluminum 


ft  6 

Modulus  of  Elasticity    29.4  x  10  psi  10.2  x  10  psi 

Yield  Point  40,800  psi  20,300  psi 

Ultimate  Strength        59,300  psi  22,250  psi 


86 


Effective  Masses 


A  very  simple  approximation  of  the  effective  weights  of  the  structure 
was  made  by  including  one  half  of  the  weight  of  two  supporting  columns 
with  the  upper  mass  and  the  weight  of  two  full  supporting  columns  with 
the  center  mass.   The  remaining  two  halves  of  the  lower  columns  were  in- 
cluded with  the  anvil  mass.   The  figure  below  shows  this  procedure. 


Ms 


M       EFFECTIVE 

i  - 


M2. 


J  M ,      F  F  FFC  T I 


M,    * 


The  values   of  the  above  effective  masses   are    listed   in  table   2. 

2 
The  units  are      lb-sec   /ft. 

Table   2 

Mass  Number  Steel  effective     Al  effective 

1  18.9  18.9 

2  0.743  0.711 

3  0.717  0.705 


APPENDIX  IV 
FURTHER  PROGRAM  MODIFICATIONS  CONSIDERED 

Two  further  modifications  to  the  computer  solution  were  considered. 
One  modification  was  the  determination  of  the  effect  of  vertical  accelera- 
tion upon  the  structure's  columns.   The  fact  that  the  anvil  experiences 
a  vertical  acceleration  upon  impact,  produces  a  force  in  addition  to  the 

structure's  static  weight.   This  vertical  acceleration  was  obtained  for 

2 
several  runs  and  found  not  to  exceed  about  5ft/ sec  .   This  would  have 

resulted  in  a  variable  correction  to  the  effective  value  of  Q  of  not 
more  than  16%,  and  the  Q/Q  correction  was  itself  rather  small,  of  the 
order  of  27».   Thus  it  was  not  consideted  of  si  f  f icient  importance  to 
make  this  correction  which  would  have  greatly  complicated  the  theoretical 
analysis  and  would  have  resulted  in  considerably  more  cumbersome  dif- 
ferential equations  to  work  with. 

The  second  modification  was  an  attempt  to  force  the  theoretical  and 
experimental  total  momentum  vs.  time  curves  to  agree  exactly,  by  varying 
the  value  of  initial  anvil  displacement  during  the  run  in  such  a  manner 
as  to  retain  the  proper  momentum  relation.   This  was  tried  for  several 
runs,  resulting  in  perfect  agreement  in  the  plots  of  total  momentum  vs. 
time;  however,  practically  no  change  in  the  velocity  response  was  noted. 
This  modification  was  likewise  not  included  since  the  constant  value  of 
a  =  0.7'  produced  essentially  the  same  results. 
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APPENDIX  V 
TYPICAL  COMPUTER  OUTPUT 
The  tabulation  of  theoretical  values  for  a  typical  run  shows  the 
time  in  MS,  velocities  V(2)  and  V(3)  in  ft/sec,  and  the  total  momentum 
in  lb-sec  designated  by  the  term  SUM, 
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Computer  Solution  STEP  s  5.02 
(T-  ^,000  psi    Steel 


J"ime         sum  vo) 

1.0000E-03                 9.481E+0J  1.312E-02 

2.0000E-03                 9.l»75fc+"6l  S.233F~-02 

?-R888E-8?                 9.U69F+01  ?!    72E-0T 

^.0000b-03                 v.4621,  +  01  2.07'OE-O 

b.0000E-03                 P.456E+01  II20-7E-0 

6.0000E-03                 y.449h*0l  *u573E-0 

7.0000E-03                 9,U43E+01  SllSOE-O 

8.0000E-03                9lu37E+0  79- 

9.0000E-03                  9.4  30L-  +  01  - 

.O000E-O2                9.U23E+01  K13UE  +  00 

•iR88I"8S  9.417E+C1  I.297E  +  00 

•1888l~2^                 ^.4 I0L  +  01  1.H53E+00 

•?888I~8?  9.403E+01  I.602E+00 

*rS22|"R?  9.3V7E+01  I.7U3E+00 

•a888,E~R?  9.390l+01  J.876E+00 

•7888^8?  g.383E*501  2.00OE+00 

1  •  7000E— 02  9.376E+01  2.11'iE  +  00 

1.8000E-02  9.369E+0I  il220F*00 

J*8888f~8?  MSI^S1  2^r6Eloo 

o*9888t"8?  v.3SSh+01  2.402E+00 

2.1000E-02  9..3U8E+01  2.479F+00 

2-?000E-02  9.341E+01  21547^00 

?'i5888f~R2>  9.3  3^L+0I  2.606E  +  00 

?"S888I~8§  g.327E+01  2.657E+00 

i*58881"8?  9-320E+01  2.699E+00 

^•£222E~02  9.313F+01  2.73&E+00 

H888iE~8§  ?.506t+01  2I763E+00 

9'o888t"8^  X-?X9£+8  2.785E+00 

^•n888rE~8?  V.291EJ  +  01  2.801E+00 

?'?88nl"8^  g.28UE  +  0J  2.8  131:  + 00 

^.1000E-02  /.276F  +  0  2    HP^F+OO 

3.2000E-02  9.269E  +  0  yltiiftoQ 

3. 3000E-02  9.26  1 L  +  01  2    83    E  +  00 

3.4000E-02  9.254F+01  7    83fcF+nn 

1-2IW8M}?  H!^*8  2    837^00 

3.6000E-02  y.239t+0  2.841E+00 

4-£222F"02  ?.231E+01  2I8U7E+00 

*-90nnP~8?  H2?£+8!  2.856EtoO 

h*n888rE"P?  9.2 16E  +  01  2-869E+00 

4.0000E-02  9.203I-  +  01  2-88?F+On 

4.1000E-02  9.200E+0:  2I910E+00 

4.2000E-02  9.193K+01  u    9 

4.3000E-02  9!    85E+,0  2.978E+00 

£-§888l"8S  V\lvA\  3loi9Etoo 

u'Annnp"8^  ?•    £?l+8  I  5.07IE  +  00 

/•2x9^-(j2  9.16    t:  +  0  3.M1M00 

u-aoon^"^  H?ff+8!     •  |-200Et8o 

m-  .  o  0  U  0  E  -  0  2  / .  1  4  t>  t-  +■  0  4    >  7  rt  r  »  n  n 

S-o888l"R§  HUH  3:366Et00 

b. OOOOE-02  '.129t+0  3.4621-fOO 

H888f~8§  H?ll+8!  SIsSIeJsS 

D.^OOOE-02  9.1  3-+0  3  An?e  +  nn 

5.3000E-02  9.  10SF+0  !  '"n    S 

S.if000E-02  9.  096b +  0  3!?  I  I"  () 

5.5000E-02  9.088F+0  J*,  , 

5.6000E-02  9  080F+G  ^228^00 

5.T000E-02  9.072E+01  ul38?E+00 

^•S822E-°2  9.063E+01  Si-  „<  too 

5.9000E-02  9.05SE  +  QH  4      +q 

6.0000E-02  V.0U7F  +  01  Sl^QE  +  OO 

! 

0 

90 


VI  3) 

5 

.  986 I 

—  06 

9 

.56  If 

-O'i 

it 

,8271 

-0*1 

1 

.'.  191. 

-0  3 

3 

.69  H 

-0  S 

( 

•  608t 

-0  1 

1 

.  W71 

2 

.  5681 

-0  2 

3 

.  r>si 

-02 

'> 

.6M9E 

-0? 

H 

.  1  3  3* 

-0,' 

I 

.  1281= 

-0  I 

1 

.5  16b 

-0| 

1 

.983F. 

-of 

2 

.S35R 

-111 

3 

.  I76F 

-0  I 

0 

.9081: 

-0  I 

'1 

.7351 

-0  I 

f> 

.6'.  7F 

-0  I 

6 

.6751 

-0  I 

f 

.  '880 

-0  I 

8 

.99  31- 

-0  I 

.02  9fc 

1  00 

1 

.  16/1. 

too 

1 

.313 1 

<(j() 

1 

.'i/,  7F. 

»()() 

1 

.  6  y  B 1 

1  00 

1 

.  fr,\-  k)o 

I 

.96  BE 

l  00 

2, 

.  IU5F  +  00 

2 

. 52SF+00 

2. 

.S09F+00 

2, 

.  6v 'n 

K)0 

2. 

,880E 

» 00 

5, 

.0661 

1  00 

3. 

.  2S  or- 

+  0(1 

5. 

.  '1  3  2  1 

1  00 

'>, 

.61  it-  ♦  00 

$< 

,  /M6I 

♦  00 

3. 

,9'»6l- 

t  00 

Ill, 

,  12  01- 

t  00 

'1. 

27  n 

♦  00 

'1  . 

u  >  n 

♦  00 

it. 

569E- 

«  00 

'1. 

70  5E 

♦  00 

'1. 

H  '  71  too 

4. 

7»t  5E  ♦  00 

'>. 

0U9I  - 

too 

9. 

l'i5b- 

>  00 

'>. 

2321  1 

1  00 

'>. 

30MI-  ♦()() 

51. 

3/'6l-»  00 

5, 

1*  -1  'l  b  i 

■  go 

'>'. 

'tM3b'l 

00 

s. 

')2'H   l 

00 

'v. 

V>6|   I 

00 

'». 

S82I   1 

00 

'■. 

60  ll  *  00 

'>. 

6  Mil  *  00 

'). 

6221  i  oij 

TIME 

6.  ioooe- o: 

6.20006- 
6. 5000E-02 
6.^0006-02 

6.50006-02 
6.60006-02 
6.7000E-02 
6.80006-02 
6.9000E-02 
7.00006-02 
7.  10006-02 
7.20006-02 
f. 30006-02 
7.40006-02 


SUM 

K   ■ 

>. 

■ .  0  2  1 

.  0  1   3 

•/.OOu 
8.995 
:•, .  9  8  7 
8.978 
8.969 
n.^'>  ; 
8.95  I 
3.94  5 
8.  934 
t.925 


F  *     ' 

H  4  0  1 

6  +  01 
6  +  01 
t-  +  01 
E  +  01 

r-+oi 

E  +  01 

e+oi 

E  +  01 
E+01 
E  t-01 
F  +  01 


V(2) 

E  +  0< 

'5 .  2  / 
..     886  +  C 
5-55  26+00 
5. 7096+00 

5.  8596+00 

5  .  ■.'  ^9  6+00 
6.  129E  +  00 
6.2466+00 

6. 5506+00 
6.4396+00 

6.3  1  36  +  00 
6.5706+00 
6.61  1E  +  00 


V(3) 

5.6276+00 
s. 6286+00 

5  .  6  '-1  /  E  +  0  0 
5.6246+00 
3.6?  1E+00 
5  .  6  1 9  E  +  0  0 
5.6 I8E+00 
5.619E+00 
5.62  36+00 
3.6  3GF+00 

5.6-;  le+oo 

5. 6566+00 
5.6776+00 

3. 7026+00 


7.5000E-02 

7.60006-02 

7.70006- 

7.80006-02 

7.9000F-C2 

8.0000E-02 

8.  10006-02 
8.20006-02 
8.3000E-02 
8.4-0006-02 
8.5000E-02 
8. 60006-02 
8.70006-02 
8.80006-02 
8.90006-02 

9.  00006-02 
7.  I000E-02 
9.20006-02 
9.30006-02 
9.  14  0006-02 
9. 50006-02 
9.60006-02 
9.70006-02 
9.80006-02 
9.90006-02 
1  .00006-0  1 


8.916E  +  01J 

8. 9066+01 

8.8976+01 

.8886+01 

! .  8  7  9  E  +  0 1 

5.8696+01 

8.8606+01 

8.  ': 

.  ;u  lt+Oi 

B. 8326+01 
...  52  >E  +  01 
8.8 126  +  I 

.  '  0  ••  I  i  I 

1.7836  +  01 
8.7  736+01 
8.7< 

8. 7546+01 
8.744 
■■■ .  7 

■■ .  j? '■■•■ 

U.  7  146  +  01 
8.7036+01 

8 .  6rt  st  +01 
.  5  736  +  01 


6.6346+00 

6.6406+00 
6.6296+00 
6. 60 16+00 
6.5366+00 
6.4966+00 
6.4-216+00 
6.  j  126+00 
6 .2316+ 
6.  !  196  +  00 
5.9976+00 
5.8666+( 
5.7296+1 

5.5886+00 

5.4U3E+00 
5.2976+00 

5.  1'>2iH 
5.009  6+ 
4.1     i-OO 
4. f366+00 
4.6096+00 
4.4  ?06+00 
4.38 16+00 
4 . 2826+00 
4.  1  94 E * 
4. 1 19  6+00 


3.7S26+00 
5.7676+00 
5.8066+00 
5.85 16+00 
5.8996+00 
5.96  16  +  00 
6.0H66  +  0  0 
6 .0636+00 
6.1215+ 
6.  1806  +  00 
6.2396+00 
6.29 re+oo 

6. 5526+00 
6.4046+00 

6.45  16+00 
6.4  ''56  +  00 
6. 3 2 6 E+ 00 
6.53  66+00 
6  .  3  7  5  E  +  0  0 
6 .  5  3  6  F  +  0  0 
6.5856+00 
6.574  6+30 
6 .  55  1 E  +  0 0 
6.3  I7E +00 
6.47  16+00 
6.4 126+00 
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Computer  Solution  STEP  =    5*. 02 
CT=   50,000  psi         Steel 


TIME 
1.000 
2.000 
3.000 
4.  000 
5.000 
6.000 
7.000 
8.000 
9.000 
1 


000 
100 
200 
300 
400 
500 
600 
700 
800 
900 
000 
100 
2.200 
2.300 
2.400 
2„500 
2.600 
2.700 
2.800 
2.900 
3.000 

3.  100 
3.200 
3.300 
3.400 
3.500 
3.600 
3.700 
3.800 
3.900 
4.000 

4.  100 
4.200 
4.300 


0E-03 
0E-03 
OE-03 
OF- 03 
OF- 03 
OE-03 
OE-03 
OE-03 
OE-03 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OF-02 
OE-02 
OF-0? 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OF-02 
OE-02 
OE-02 
OE-02 
OF-02 
OF-02 
OE-02 
OE-02 
OF-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 


4.400 
4.500 
4.600 
4.700 
4.800 
4.900 
5.000 
5.  100 
5.200 
5.300 
5.400 
5.500 
5.600 
5.700 
5.800 
5.900 
6.000 


OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 
OE-02 


SUM 

9 . 4  8  U  +  0 
9. 4 75 E+0 
/.467E+0 
9.462E+0 
9.4  56F+0 
V. 44V 1+0 
/  .  4  4  S  r  +  0 
9.43/E+O 
9.430E+0 
>.423h+0 
y .  4  1  7  E  +  0 
v.'i  10E  +  0 
9.4  04F+0 
7  .  3  9  ,'  E  +  0 
J.390E+0 
9.383E+0 
J.377E+.0 
>. 370E+0 
J.363Q+0 
9.3566+0 
9.349E+0 
9.342E+0 
>.335E+0 
9.328E+0 
9.  32  IE  +H) 
'.  S  14 E+0 
7.  S06t+-0 
7.2 97 F+O 
7.292E+0 
9.285E+0 
9 .  2  77  E  +  0 
/.27GE+.0 
9.262E-M0 
9.255E*0 
9.2  48F+0 
9.240E-HO 
V.233l+C 
7.225F+C 
c/.2  17t-t-0 
9.2  10E  +  0 
9.202E+.C 
9.  194 E+0 
7.  1 8  7  r  +0 


179 
171 
163 
155 
1  U  7 
137 
131 
123 
1  15 
10/ 


9.099 


091 
083 
0  75 
066 
0  58 
0  50 


E  +  0 
£  +  0 
t  +  0 
E+0 

L  +  0 

E+d 

E+0 
l:+'0 
E  +  0 
F  +  O 
E  +  0 
F  +  O 
E  +  0 
I  +,0 
F+'O 
E  +  0 
E  +  0 


V(2) 

1.312 
5.2  33 
1  .  1  (? 
2.010 
3.20/ 
4.574 
6.150 
7.724 
9.876 
1.  198 
1  .4  OH 
1.611 
1  .805 
1  .79  I 
2.167 
2.333 
2.488 
2.651 
T6U 
88-5 
794 
072 
I  79 
254 
3  210 
3.375 
3.42? 
3.1-160 

3.4  70 

3.5  14 
3.S  S3 
3.54  6 
3.5  57 

3.5  65 
3.5/2 
3  .  5  n 
3.58H 
3.')  9  7 

3.6  14 
3.  6  S3 
3.658 
3.690 
3.  730 


E-02 
E-02 
F-0  I 
E-0  1 

I; -01 

F-0  1 
E-0  1 

E-0  I 
E-0  1 
E  +  00 
F+00 
E  +  00 
E  +  00 
E  +  00 
E  +  0  0 
E  +  00 
E  +  00 
E+00 
F  +  00 
£♦00 
E+00 
E+00 
E  +  00 
£♦00 
F  +  00 
E+00 
E  +  00 
£♦00 
£♦00 
£♦  0  0 
£  +  00 
I:  ♦  0  0 
£  ♦  0  0 
E  +  00 
E  +  00 
E  +  00 
£♦00 
£  +  00 
£  +  00 
E+00 
E  +  00 
E+00 
E  +  00 


3.778E 


7  05E 
9  HOE 

0  69E 

1  6  9  £ 
280  E 
40  3 £ 
537F 
68I2E 

8  S|6  £ 
97  8  E 

5.  167E 
5.34'1E 
5.51  8 1 

5.67/1 
5.0/6E 


+  00 
+  00 
+  00 
+  00 
+  00 
+  00 
+  00 

♦  00 

♦  00 
+  00 
+  00 

♦  00 
»  0  0 
«  00 

♦  00 
+  00 

♦  00 


.9861 

56  II 
M2  f\ 

5  191 

6  7  1 1 
608 

3  7-91 
S  6  H I 
7571 
6671 
20  1 1 
1461 
5541 

0'>7l 
64  f\ 

1471 
1561 
Offll 

1  161 
?72l 
5461 
9  WE 

144 
50  61 

4  781 

66  11- 

R52H 

05  SI 
26  II 
4761 

6761- 

7  7  11 
1501 
SMOI 

6  |  2  £ 
H»l  51 
0  73t 
30  IE 
5  7  4  F 
74  _S£ 
9  5  5  £ 
16  IE 
358E 

5.547E 
5./2  5E 
H9UE 
0 5  I  F 
I98E 
SVE 
4V>£ 

5  6  6  F 

6  6  5E 
f52E 
82  8E 
FJ  }  SI. 
7 'i  8E 

7  7  S| 
0  7  7 1- 
05  (\ 
078E 


2 


6. 
/. 
8. 
9. 

1  . 

1. 

1. 

I. 

1. 

7. 

7. 

7. 

2. 

2. 

3. 
j 

>  • 

> . 
S. 
4. 
4. 
4. 
4.. 
4. 
5. 
5. 


-  0 


-  (I 


-0 
1-0 


I  1 1 


-0 
-  0 


-0 
-0 
-0 

~l 

-  t 

I 

i  ( 
« ( 

♦  I 
i  t 

♦  I 
t  I 
1 1 
» I 
1 1 

I  ( 
1 t 

I I 
» I 
+  ( 
+  0 

+  ( 

+  0 
+G 

+  ( 
+  0 
+  0 

+  1 
+  ( 

1 1 

♦  I 

I  ( 

I I 
• 
* 
( 


(I 
!. 
( 

+  0 

+  ( 
» ( 

♦  ( 

♦  ( 

» I 
»  I 
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tt:  e 

sT: 

T-r(2) 

6.  I000F-G2 

) .  0  4  1     ■ 

6. 

6.2000E- 

■ .     :    . 

■     6. 

6.3000E-0; 

6 .     • 

6.4000F-02 

■>.0  lo 

6. 

6.5000E-02 

9.007E+0 

6- 

6.6000E-02 

- . 999E+0 

6.844 

6. rOOOF-02 

.  )  -I  \      i    | 

6  .971 

6.8000E-02 

7. 

6.9000E-02 

i .  9  7  3E  +  0 

7.  184 

7.0000F-02 

.  ?64F+0 

7.266 

7.  1000E-U2 

7.3 

7.2000E-02 

8.946E+0 

7. 

7.3000E-02 

a.  9  3  re+o 

7. 

7.4000E-02 

.928E+0 

7  .  u  l  8 

7.5000E-02 

8.919E+Q 

7 . 

7.60006-02 

.910E+0 

7 . 

7. 70006-02 

8.  '•()  IE  +  Q 

7.335 

7.8000F-02 

- .891E+0 

1                  7.270 

7.9000E-02 

. 832E+0 

r. 

8.0000E-02 

8.873E+0 

7.1 

8. 1000F-02 

.8631  +0 

6.976 

«-2000E-02 

-  8  54 E  +  0 

6 .  bU  8 

8.30006-02 

. 8446+0 

6. 

8.4000E-02 

.835E+0 

6.556 

8.5000E-02 

!  .8 2  i  E + 0 

6.394 

8.6000F-C2 

.  3  1  5  -  +  n 

6.  224 

8.7000E-02 

8. 805 t+0 

6.04  8 

8.8000E-02 

i.  f96E+0 

5. 

8.9000E-02 

.  7  86»-  +  0 

1                 5. 

9.0000E-02 

i.7766  +  0 

5.1 

9.  IOOOE- ;.  i2 

.  7661  + 

5. 

9.2000F-U2 

.  f !  i  < 

1                 5.139 

9.3000E-02 

8.  rt+6i  +0 

1                 4.964 

9.40006-02 

.736E+0 

4. 

9.5000E-02 

8  .725 

4.< 

9.6000E-02 

.  n  5  E + 0 

I                  4. 

9.7000E-02 

• .  rose +  : 

1                    4  . 

V. 800  0 E -02 

I.695E  +  0 

1 

9.9000E-02 

i.68*i     • 

1                 4.087 

1  .OOGOF-01 

.  :>  (  4  ■:  +  o 

1                  5.980 

E  ^  0  0 
E+  0  0 
b  +  00 
Ft- 00 
E  +  OC 
E+00 
E  +  00 
6  +  00 

E  +  00 
E  +  00 
E  +  00 
E+00 
E  +  00 
F  +  00 
E  +  00 
F  +  00 
E  +  00 
E  +  00 
E  +  00 
6  +  00 
E  +  00 
E  +  00 
F  +  00 
E  +  00 

6  +  00 
6  +  00 
6+00 
6+00 
6  +  00 
E  +  00 
E  +  00 
E  +  00 
E  +  00 
6  +  00 
6  +  00 
£+00 


7(3) 

0  9  9 


102 
0  ;v 
093 
r.086 
7.0 
f.0  7G 
7.C 
7.C 
7.05  3 
7.050 
7.0 'I  9 
7.0  GO 
7.053 
7.1 


065 
0  72 

I 

0  •'  7 
103 

10  7 
108 

0  )  8 

064 

OS/ 

0  0  1 
955 

f5  8 

56  9 
4  5  6 

3  ;>  i 

1  9  4 


6  +  0  0 
t-  +  00 
6  +  00 
E  +  OC 
F+00 

E  +  0  0 
6  +  00 
6  +  00 
F  +  00 
E+00 
E+00 
E  +  00 
6  +  00 
F  +  GO 
6  +  0  0 
6  +  00 
F  +  00 
6  +  00 
E  +  00 
6  +  00 
E  +  00 
E+00 
F  +  0  0 
F  +  00 
6  +  00 
6  +  00 
E  +  00 
E  +  00 
6  +  :)0 
6  +  00 
6  +  0  0 
\}  +  0V 
F  +  00 
F  +  00 
E  +  00 
6+00 
F  +  00 
6  +  0  0 
6  +  0  0 
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Computer  Solution  STEP  =  5.02 
<T=  60,000  psi    Steel 


TIME 

SUM 

VI?) 

V{  3  ) 

1.00006-03 

9.48 1E+0 

1  .3  121.-0? 

5.  9  8  A  6-  Oft 

2.00006-03 

7.47SE+01 

5. 23  3  6-02 

9.36  lt-l)S 

3.0000E-03 

9.4696+01 

1.  172E-01 

4  .82  (t-n>\ 

4.0000E-03 

9.462E+01 

2.070E-0 1 

I.519E-0  '. 

5.0000F-03 

9.456E+0 

3.2076-01 

"..6/U  -0  3 

6.000GE-03 

9.4496+0 

4.573E-0  1 

7;  6086- 01 

7.0000E-03 

9.443E+01 

6.  I50E-0 1 

1  .  3991  -0? 

6.0000E-03 

9.4  37E+0 

1                      7.924E-01 

2.  3681  -0? 

9.0000E-03 

9.43OE+0 

9.8766-01     • 

i./sn  -0? 

1.0000E-02 

9.423E+0 

1                          1  .  1981-  +00 

5.66  f\  -0? 

1  .  1000E-02 

9 .  4  1  7  E  +  0 

1                      1.423E+00 

8.2026-0? 

1.2000E-02 

9.4  10E+0 

1                       l.oSBE+00 
1 .6996-1  00 

1.  Pi  7  6-0  1 

1.3000E-02 

9. 4 04 E +01 

1.5591  -0  1 

1  .4000E-02 

9.397E+01 

2.  1  326+00 

2.0651 -0 1 

1. 50006-02 

9.3906+01 

2.355E+00 

2.6^7E-0  1 

1.6000E-02 

9. 3846+01 

2.566E+00 

3  .  4  0  ',  I  -  1 )  1 

1.7000E-02 

I    9.377E+01 

j                  2.765£+00 

'1  .?'>  IE— 0  1 

1.8000E-02 

9.370E+01 

2.9521E+00 

5.2271  -0  1 

1.9000E-02 

9. 36'rii  +01 

3.  125E+00 

/>.  3  36  6-0  1 

2.0000E-02 

V.356F+01 

3.2866+00 

r.580E-0  1 

2.  1000E-02 

9.349E+0 

3.43  36+00 

8.96 iE-0 1 

2.2000E-02 

9. 342F+01 

3.566E+00 

1  .1)496-1  00 

2.3000E-02 

9.335E+0 

3.686E+00 

1  -  ?  Pi  6*00 

2.4000E-02 

9.328E+0' 

3.79  36+00 

1  .  394  6+00 

2.50006-02 

9.3216+0 

3.88  76+00 

1  LS86I  1  00 

2.6000E-02 

9 . 3 1 4  E+ 0 

3.9696+00 

1  '.  792F+O0 

2.7000E-02 

9.307F+0 

4.04QE+00 

4.  1006+00 

2  S.009F+00 

2.8000E-02 

9.300E+0 

?;237F+00 

2-9000E-02 

9.293E+0 

4.  H+9E  +  00 

2l4  736*  00 

3.0000E-02 

9. 2 856 +0 

4.  1906+00 
4.2236+00 

>»'.  ( .>  31   100 

3.  1000E-02 

9. 2 7 BE +0 

2 '.9786  1  00 

3.2000E-G2 

9.27 1E+0 

4.250E+O0 

3.2  4  IE +  00 

3.3000E-02 

9.263  6+0 

!                       4. 27  IF  +  00 

3.3081   1  00 

3.4000E-02 

7.2  566+0 

4.2.88E  +  00 

3. 7801  *  00 

3.5000E-Q2 

9.249E+0 

1                      4. 3036+00 

4.0331    »00 

3.6000E-02 

9.2416+0 

4.3  166+00 

4.33  1!  +  00 

3..7000F-02 

9.234E+0 

1                       4.3296+00 

4.6086* 00 

3.8000E-02 

9.226E+0" 

1                       4.3446+00 

M.88?6+00 

3.9000E-02 

9.219E+0 

4. 36  J E  +  00 

S.  IS 'if  +00 

4.0000E-G2 

9  .  2  1  1  E  +  0 

1                       4.38  36+00 

5.U2  11+00 

4. 1000F-02 

I     9.204E+0 

1                       4.U 106+00 

5.6831 +00 

4.2000E-02 

9.  1966+0 

4.4446+00 

5.9  371  +00 

4..3000E-02 

9.  I8H1-  +  O 

4.485E+00 

6.  1846  +  00 

4..4000E-02 

9  .  1  8  1  E  +  0 

1                        4.5356+00 

■      6.42  If  +00 

4.5000E-02 

9.  173E  +  0 

1                        4.S9SL+00 

6.6486-1  00 

4.6000E-02 

'.  165E+0 

I                       4.665E+00 

6.86  31  +o() 

4.7000E-02 

9.  157r+0 

1                        4. 744 Fi 00 

7.06  71  *oo 

4.8000E-02 

9.  1496+0 

i                         4.8346*00 

7.?3Bi  +00 

,4.9000E-02 

9  .  1 4  2  E + 0 

1                         4.9336*00 

/.'li/.i  « 00 

5.0000E-02 

9. 1 346  +  0 

1                         5. 0466  +  00 

7.  6i)OI+00 

,5.  10006-02 

9.126-r  +  O 

1                         5.  1666  +  00 

7.  7'>0I  +00 

.5.20006-02 

9.  1  18F  +  0 

1                       5.2966+00 

7.8MM  +  00 

S.3000F-02 

9.  1  IOl+0 

1                       5.4346+00 

8*,00/E  +  00 

5.4000F-02 

■>.  I02E+0 

5.58QF+00 

Hi  1  1ST  ♦•()() 

5.5000F-02 

9.094F+0 

t                       5. 7336+00 

8w209E+00 

5.6000E-02 

9.085E+0 

1                      5.89JE+00 

8i.289F+O0 

5.700OE-02 

9.0  771; +0 

1                       6.05*6*00 

8.  -I'./.l  +00 

5.8000E-02 

9.0  69E+0 

1                       6.? IMF +00 

M.4  1  IF  +00 

5.9000E-02 

9 .  0  6  1 1  +  0 

i                      6.3B4E+00 

H.llS  31+00 

6.0000E-U2 

9.053E+0 

1                       6.SS0E+0O 

H.4BSF+00 
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TIME 
6. 1000E-02 

6.2000L- 

6 .  3 

6.4000E-02 
6.5000F-G2 
6.6000F-02 
6.7000E-02 
6.8000E-02 
6..9000L--02 
7.0000L-02 

7.  1000E-G2 
7.2000E-02 
7..3000E-02 
7. 4000 E- 02 
7.5000E-02 
7.6000E-02 
7.7000E-02 
7.8000F-02 
7.9000E-02 
6.0000E-02 

8.  1000F-02 
8.200UF-02 
8.30OQE-02 
8.4000E-02 
8.5000E-02 


SUM 

V(2) 

V(3) 

9  .  0  4  U  !■  f  0 1 

6  .  7  1  3  E  ■+  0  0 

■ .505E+00 

>.0  U  1  «-01 

.    6.8731 

-  .5  I6E+00 

9.027!  +  (  " 

7.0  !8E+i  . 

8.5  18E-t 

9.0 J9E+01 

7.  17 'j  F +  00 

8.512E+00 

9  .  0  1 0  E  +  0 1 

7.314E+00 

8.498E+00 

9.002F+     ■ 

7.443E+00 

.478E+00 

■  .  •  -  5    '  01 

7. 560 E+00 

8.452E+00 

8.984E+01 

7.663E+00 

8.422E+ 

I.976E-I  01 

7. 753E+00 

8.38  7E+00 

^.9671  +01 

7.826E+00 

3.348E+00 

.  •)"■  hOl 

7.^;;3F  + 

8.307E+ 

8.949E+01 

7 . 9  2  i  E + 0  0 

8.264E+00 

• . 940F+01 

7* 944 E+00 

.219E+00 

.  .931E+01 

7  . 9  4  7  E  +  0  0 

3 . 1 7  2  E  +  0 0 

8.922F+01 

7.931 E  + 

. I24E+00 

8  .  9  1 3  E  +  0 1 

7.895E+C 

■  .074E+00 

3.903E+01 

7.84 1E+00 

8. 024 E+00 

1.894E+01 

7.769E+00   . 

7.973E+00 

.      15    +01 

7.678E+00 

7  . 9  ?  1 E  + 

H.875E+01 

7.570E+00 

7  .  8  h  7  F  -i 

8.866E+01 

7.44 5 E+ CO 

7.812E+ 

8.856E+Q1 

7. 505E+00- 

7.755E+00 

8.847E+01 

7.  150 E+00 

7.695E+00 

8.837E+01    • 

6.982E+0 

7.632E+00 

8.827F+01 

6.^03E+00 

7.56  6F+00 

8.600 
8.700 
8.800 
8.900 
9.00  0 
100 
200 
300 
1+00 
9.500 
9.600 
9.700 
9.800 
9.90  0 
1  .000 


OE-02 

OF-02 
OE-02 
OE-02 
-02 
OF-02 
OE-02 
I  I  -02 
31  -02 
'  I  -02 
0E- 
OE-02 
OE-02 
0|  -02 
OE-0  1 


e 

.  8  1  7  E  +  0 1 

a 

.808E+01 

.798E+01 

•' 

.788E+01 

■ 

.77RE+01 

• 

•768E+01 

.  758E  +  01 

H 

.7471  +01 

8, 

.73  71-  + 01 

.' 

.72  7F  f-OI 

.  7J7E  +  01 

.  7061  fOl 

.696E+01 

•  686F  i-OI 

• 

.  6  75 1 +  01 

6 .  6  1  4 
6.4  16 
6.212 
6.012 

5.  ne 

5  .  6  3  2 

5. 

5.282 

5.11V 

M.  i>65 
■•  .820 

'i  .  6  8  3 
4.554 

U.U  '!i 

4.  j?2 


F  +  00 
E+00 
E+00 
fc  +  OO 
E  +  00 
E  +  00 
E  +  00 
E  +  00 
t  +  00 
E+00 
E+00 
E+00 
E  +  00 
E+00 
E+00 


7. 
7. 

r, 

7 

7! 

7, 

6. 

6. 

6. 

6. 

6 

6 

6 

5 

5 


49 
42 

33 
25 
15 

05 

9  4 
8  3 

57 
42 
2  7 
1  1 
9U 
7  7 


5  E  +  00 
OE  +  00 
9  E  +  0  0 

2E  + 
9  E  +  0  0 
8E  +  00 

Q  F  +  0  0 
2E+00 

6F  +  00 
2E  +  00 
9  E  +  0  0 
7  F  +  0  0 
7  F  +  00 
8E  +  00 
OE  +  00 
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APPENDIX  VI 
INSTRUMENTATION  SERIAL  AND  MODEL  NUMBERS 


Component 

Model 

Velocity  meter  on  M1 

* 

Velocity  meter  on  M» 

* 

Velocity  meter  on  M„ 

* 

Accelerometer  on  M~ 

A5A-3OO-50O 

Accelerometer  on  M 

A5A-300-500 

Velocity  Meter  Control 

122C 

Unit 

Oscillograph 

5-119 

Amplifier  #1 

1-113B 

Amplifier  #2 

1-113B 

Power  Supply 

2-105 

1000  cycle 

723-C 

Serial 
159 
222 
109 
4369 
4012 
101 

109025 
198BL10 
641BD09 
245  R4 
2326 


*None  available 
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